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Abstract.  Theoretical  analysis  of  the  spatial  hole  burning  in  quantum  dot  (QD)  lasers  is  given. 
The  multi-mode  generation  threshold  is  calculated.  The  processes  of  the  thermally  excited  escapes 
of  carriers  away  from  QDs  are  shown  to  control  the  multi-mode  generation  threshold.  The 
dependences  of  the  multi-mode  generation  threshold  on  the  root  mean  square  of  relative  QD  size 
fluctuations,  cavity  length,  surface  density  of  QDs,  and  temperature  are  obtained. 


Introduction 

In  [ 1]— [3],  theory  of  threshold  current  density  of  a quantum  dot  (QD)  laser  and  its 
temperature  dependence  has  been  developed  having  regard  to  inhomogeneous  line 
broadening  caused  by  the  dispersion  in  QD  sizes.  The  optimum  parameters  of  the  laser 
structure  minimizing  the  threshold  current  density  have  been  calculated  as  the  functions 
of  the  QD  size  dispersion,  total  losses,  and  temperature. 

This  article  discusses  the  effect  of  spatial  hole  burning  and  multi-mode  generation 
threshold  in  QD  lasers.  As  in  conventional  quantum  well  (QW)  or  bulk  lasers  (as  well 
as  in  solid  state  lasers)  [4,  5],  spatial  hole  burning  in  QD  lasers  is  due  to  the  non- 
uniformity of  the  stimulated  recombination  of  carriers  along  the  longitudinal  direction 
in  the  waveguide.  Due  to  the  fact  that,  at  and  above  the  lasing  threshold,  the  electric 
field  of  the  emitted  light  is  a standing  wave  and  is  a periodic  function  of  the  longitudinal 
coordinate,  the  stimulated  recombination  of  the  carriers  will  be  more  intensive  in  the 
QDs  located  at  the  antinodes  of  the  light  intensity,  while  it  will  be  less  intensive  in  the 
QDs  located  at  the  nodes.  As  a result,  overfilling  of  the  QDs  located  near  the  nodes 
may  take  place.  This  leads  to  the  lasing  generation  of  the  other  longitudinal  modes 
(together  with  the  main  mode)  with  antinodes  distinct  from  those  of  the  first  mode.  A 
problem  of  the  multi-mode  generation  is  of  first  importance  for  the  laser  applications. 
A study  of  the  physical  processes  controlling  the  multi-mode  generation  threshold  is 
necessary  to  find  the  ways  of  suppressing  the  additional  modes  and  to  offer  the  proper 
design  of  single-mode  operating  lasers. 

1 Processes  controlling  the  spatial  distributions  of  carriers 

In  QW  or  bulk  lasers,  diffusion  in  the  active  region  will  tend  to  smooth  out  the  non- 
uniform  carrier  distributions  and  population  inversion  along  the  longitudinal  direction, 
thus  suppressing  totally  or  partly  the  effect  of  spatial  hole  burning  [4,  5], 

In  QD  lasers,  diffusion  will  play  a similar  yet  minor  role.  The  point  is  that  the  carriers, 
contributing  to  the  stimulated  emission,  are  those  totally  confined  in  QDs.  There  are  also 
free  carriers  in  the  optical  confinement  layer  (OCL)  which  contribute  to  the  spontaneous 
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emission,  thus  increasing  the  threshold  current  density.  The  free-carrier  densities  and 
the  confined  carrier  level  occupancies  in  QDs  are  coupled  to  each  other  by  the  rate 
balance  equations.  Due  to  this  coupling,  diffusion  of  free  carriers  should  equalize  to 
some  extent  the  level  occupancies  in  different  QDs.  Hence,  two  processes  control  the 
spatial  distribution  of  free  and  confined  carriers  along  the  longitudinal  direction.  These 
processes  are:  the  thermally  excited  escapes  of  the  carriers  from  QDs  to  the  continuous 
spectrum  states  and  the  diffusion  of  free  carriers  along  the  longitudinal  direction. 

The  slowest  process  of  the  two  above  controls  the  carrier  space  distribution.  In  this 
work,  the  thermally  excited  escapes  from  QDs,  rather  than  the  diffusion,  are  shown 
to  limit  smoothing-out  the  carrier  space  distribution.  There  is  an  evident  analogy  to 
spatial  hole  burning  in  bulk  lasers  containing  impurity  centers  in  the  active  region  [6,  7], 
Non-vanishing  values  of  the  characteristic  times  of  thermally  excited  escapes  are  shown 
to  control  the  multi-mode  generation  threshold  in  QD  lasers. 

The  following  simple  reasoning  is  worth  presenting  here.  The  fluxes  of  the  thermally 
excited  escapes  of  the  electrons  and  holes  from  QDs  are  proportional  to  /n.pA's/'nf.p, 
where  /n.p  are  the  mean  electron  and  hole  level  occupancies  in  QDs,  N$  is  the  surface 
density  of  QDs,  t|iP  are  the  characteristic  times  of  thermally  excited  escapes  of  electrons 
and  holes  from  QDs  being  given  as  [1] 

rg  = — — = — 1 — Ts  = — = — 1 — . (1) 

” 7n«l  crnVn«l  P 7p/>l  fpVpPi 

Here  7„iP  = crn,pVn,p,  crnjP  are  the  cross  sections  of  electron  and  hole  capture  into 
a QD,  and  vnjP  are  the  thermal  velocities  of  electrons  and  holes.  In  Eq.  (1),  n\  = 
A(?CL  exp[-(A£,c  - en)/T]  and  p\  = A^CL  exp[— (Aisv  — ev)/T]  where  = 

2(m°$LT /IrrH2)^2  are  the  conduction  and  valence  band  effective  densities  of  states 
for  the  OCL  material,  AEC  and  A £v  are  the  conduction  and  valence  band  offsets  at  the 
QD-OCL  heteroboundary,  sn  p are  the  quantized  energy  levels  of  an  electron  and  hole  in 
a mean-sized  QD  (measured  from  the  corresponding  band  edges),  and  the  temperature 
T being  measured  in  terms  of  energy. 

The  free-hole  diffusion  flux  is  proportional  to  2 kDpp,  where  p = pi/p/(l  — /p)  is 
the  free-hole  density  [1,2 ],k=  (2ir/ A0 ) y/e,  \0  is  the  wavelength  at  the  maximum  gain, 
e is  the  dielectric  constant  of  the  OCL,  and  Dp  is  the  hole  diffusion  constant.  Since  Dn 
is  greater  than  Dp,  the  free-electron  diffusion  is  not  the  limiting  factor. 

The  ratio  of  the  hole  escape  flux  to  the  diffusion  one  is  (1  — fp)(NscrpVp) / (2kDp).  It 
is  controlled  by  the  cross  section  of  hole  capture  into  a QD  and  by  the  surface  density 
of  QDs.  It  is  typically  much  less  than  unity.  What  this  means  is  the  process  of  thermally 
excited  escapes  is  the  slowest  and  hence  the  limiting  one.  Should  this  ratio  be  close  to 
or  greater  than  unity,  both  the  processes  above  will  control  smoothing-out  the  spatial 
distribution  of  carriers.  The  multi-mode  generation  threshold  in  this  case  will  be  less 
than  that  in  the  case  of  small  ratio. 

2 Multi-mode  generation  threshold 

An  examination  of  the  problem  yields  the  following  equation  for  the  excess  of  injection 
(pump)  current  density  over  the  threshold  current  density  of  the  main  (closest  to  the 
maximum  of  the  gain  spectrum)  mode  required  for  oscillating  the  next  longitudinal 
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mode: 

c._  • • _ fel eNs fel eNs (J) 

J n Jt h gmax  (1_/n)  (l_/p)  g-xTg(1_/n)  + Tg(1_  /p)  U 

7n«1  IpPl 

where  jth  and  jj  are  the  threshold  current  densities  of  the  main  and  the  next  longitudinal 
modes,  respectively,  g™1**  is  the  gain  spectrum  maximum,  and  /njP  are  the  mean  electron 
and  hole  level  occupancies  in  QDs  required  for  the  lasing  of  the  main  mode.  The 
absolute  value  of  the  difference  in  the  gain  of  the  main  and  the  next  modes  is 


m 


c?g 

dE2 


m1  = \ 


d2g 


dE 2 


c 7 r 

'jn. 


(3) 


where  the  derivative  is  taken  at  E = E0,  E0  is  photon  energy  of  the  main  mode, 
SE  = He/  y/e )(n/L ) is  the  separation  between  the  photon  energies  of  the  neighbouring 
modes  (Am  = ±1),  and  /.  is  the  cavity  length. 

For  Gaussian  distribution  of  relative  QD  size  fluctuations  [1], 


M&l(  \ 

gmax  2 1 (Ae)inhom  J 


(4) 


where  (Ae)inh0m  = (<?n£n  +?p£p)^  is  the  inhomogeneous  line  broadening  due  to  fluctu- 
ations in  QD  parameters  (e.g.,  sizes),  qn^  = — {d\ne^p/ d\na)  and  S is  the  root  mean 
square  (RMS)  of  relative  QD  parameter  (size)  fluctuations  [1]. 

The  threshold  current  density  of  the  main  mode  is  [ 1]— [3] 


eNs,,  , un  /n/p 

7th  = —/n/p  + ebBnipi  — — — 

TQD  (1  ,/n)(l  ,/p) 


(5) 


where  tqd  is  the  radiative  lifetime  in  QDs,  b is  the  OCL  thickness,  and  B is  the  radiative 
constant  for  the  OCL. 

The  relative  excess  of  injection  current  density  over  the  threshold  current  density  of 
the  main  mode  required  for  oscillating  the  next  longitudinal  mode  is 


5-l 

7 th 


fel  TQD 
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(6) 


We  shall  restrict  our  consideration  to  the  case  of  charge  neutrality  in  QDs  when 


/n=/p 


(7) 


where  Nfm  = (4/^)(Ve/Ao)2TQD  ((Ae)inh0m/^)  P (a/r)  is  the  minimum  surface  den- 
sity of  QDs  required  to  attain  lasing  at  given  losses  (3  and  inhomogeneous  line  broadening 
(Ae)inhom  [1,  2],  £ is  a numerical  constant  appearing  in  QD  size  distribution  function, 
a is  the  mean  size  of  QDs,  and  F is  the  optical  confinement  factor  in  a QD  layer  (along 
the  transverse  direction  in  the  waveguide). 
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Fig  1.  Relative  multi-mode  generation  threshold  and  the  minimum  threshold  current  density  of 
the  main  mode  versus  the  RMS  of  relative  QD  size  fluctuations  5. 


3 Results  and  discussion 


Analysis  of  (6)  shows  that,  for  the  structure  optimized  at  a given  (Ae)inhom  (for  which 
the  threshold  current  density  of  the  main  mode  is  a minimum),  8 jl  jt h oc  [(Ae)ini,om]~2 
for  small  (Ae)inhom,  and  8j/jth  oc  [(Ae)inhom]_7/4  for  large  (Ae)inhom. 

As  is  easy  to  see  from  (1)  and  (6),  8j/jlh  depends  strongly  on  the  cross  sections 
of  electron  and  hole  capture  into  a QD,  o-n,p.  Calculation  of  a-n,p  is  beyond  the  scope 
of  the  present  article.  Here,  to  estimate  the  multi-mode  generation  threshold,  we  take 
crn  p = 10-13  cm2  (which  is  much  less  than  the  geometrical  cross  section  of  a QD).  At 
room  temperature  and  at  L = 500  pm,  Fig.  1 shows  the  relative  multi-mode  generation 
threshold  (the  solid  curve)  and  the  minimum  threshold  current  density  of  the  main 
mode  (the  dashed  curve)  versus  the  RMS  of  relative  QD  size  fluctuations  5.  Each  point 
on  the  curves  corresponds  to  the  specific  structure  optimized  at  the  given  5.  For  the 
structures  with  8 = 0.05  and  0.1,  8j / _/th  « 8%  and  2%,  respectively;  the  minimum 
threshold  current  density  values  are  14  and  25  A/cm2,  respectively.  The  inclusion  of 
violation  of  the  charge  neutrality  in  QDs  [2]  will  enhance  the  multi-mode  generation 
threshold. 

Thus  an  increase  in  the  QD  size  dispersion  not  only  increases  the  threshold  current 
density  but  decreases  the  multi-mode  generation  threshold  as  well. 
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